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USAARDC  Ballistic  Resea "  h  Laboratories 
Aberdeen  Proving  Ground,  ^ilryland 

IS.  ABSTRACT  “ 

■  A  gas  in  a  tube  is  excited  by  a  reciprocating  piston  operating  at  or  near 
a  resonant  frequency.  Damping  is  introduced  into  the  system  by  three  means t 
radiation  of  energy  from  one  end  of  the  tube,  rate  dependence  of  the  gas,  and 
boundary  layer  friction.  A  lumped  damping  coefficient  is  defined  to  account 
for  the  effects  of  all  three.  It  is  shown  that  in  the  small  rate  limit  the  signal 
in  the  periodic  state  suffers  negligible  distortion  in  one  travel  time,  and 
hence  its  propagation  according  to  acoustic  theory  is  valid.  The  shape  of  the 
signal  is  determined  by  a  nonlinear  ordinary  differential  equation.  The  small 
rate  condition  provides  a  test  of  the  applicability  of  the  theory  to  given  experi¬ 
mental  conditions.  When  there  is  no  d^ping^  shocks  are  a  feature  of  the  flow 
for  frequencies  in  the  resonant  band.  For  a  given  amount  of  damping  an  upper 
bound  on  the  piston  acceleration  is  given  to  ensure  a  shockless  motion.^  The 
resonant  band  is  analyzed  for  both  damped  and  inviscid  cases.  The  predictions 
of  the  theory  are  compared  with  experiment. 
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1 .  Introduction 

In  this  paper  we  discuss  the  periodic  vibrations  which 
result  when  a  column  of  gas  in  a  Kundt's  tube  is  driven  by  a 
piston  oscillating  at  a  near  resonant  frequency.  The  basic 
experimental  observations  are  well  documented,  (Saenger  and 
Hudson  (1960)).  When  the  piston  frequency  is  in  a  band  about 
a  resonant  frequency,  the  amplitude  of  the  response  is  markedly 
higher  than  the  piston  amplitude  and  shock  waves  appear  in  the 
flow.  These  phenomena  have  been  extensively  investigated  in 
the  recent  literature,  (Betchov  (1958),  Saenger  and  Hudson  (1960), 
Chu  and  Ying  (1963),  Chester  (1964),  Mortell  (1971a, b),  Collins 
(1971)).  Nevertheless,  there  are  associated  phenomena  rfhich 
need  further  investigation,  and  some  aspects  of  the  various 
analyses  which  need  clarification.  Two  questions  which  are 
resolved  concern  the  range  of  validity  of  the  usual  modifica¬ 
tions  of  acoustic  theory  used  previously,  and  the  existence  of 
a  critical  amount  of  damping  which  ensures  a  shockless  motion. 

The  basis  for  the  analysis  given  here  is  the  fact  that  the 
motion  of  the  gas  can  be  represented,  to  first  order  in  the 
amplitudes,  as  the  superposition  of  two  noninteracting  simple 
waves  traveling  in  opposite  directions,  (see  Mortell  and  Varley 
^.1970)).  This  implies  that  the  travel  time  of  a  component  wave 
in  the  tube  is  determined  by  its  own  amplitude,  and  then  any 
distortion  of  a  signal  is  self-induced.  The  essence  of  this 
approximation,  which  is  the  first  term  in  a  regular  perturbation 
expansion,  is  that  while  the  amplitudes  remain  small,  ju]  <•<  a^  , 


f 
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I  the  acceleration  in  the  gas,  ,  is  not  restricted.  This 

contrasts  with  linear  acoustic  theory  where  both  juj  << 
and  If  <<  are  required  for  its  validity,  (Here,  u  is 
the  particle  velocity,  a^  is  the  equilibrium  sound  speed  and 
L  is  a  typical  length  of  the  medium.)  A  nonlinear  theory  such 
as  this  has  been  used  for  many  years  in  gas  dynamics  when  the 
disturbance  is  generated  by  the  passage  of  a  single  progressing 
wave  (see  Whitham  (1952)).  The  noninterscting  simple  wave 
representation  ran  be  viewed  as  a  generalization  o-p  ’.‘.ftitham's 
rule  (1952)  to  wave  motions  having  two  components. 

We  consider  the  motion  of  a  gas  in  a  tube  which  is  driven 
by  a  piston  oscillating  at  one  end.  The  other  end  is  ''partially 
open"  in  the  ssnse  that  the  system  is  considered  to  lose  energy 
by  radiation  through  this  end  into  the  adjacent  medium.  We  seek 
the  time -periodic  response  of  the  gas  to  these  boundary  condi¬ 
tions.  Using  the  simple  wave  representation  the  problem  of 
calculating  the  shape  of  the  signal  at  a  boundary  is  reduced  to 
finding  solutions  of  a  nonlinear  functional  difference  equation. 
The  signal  may  distort  as  it  travels,  but  its  shape  at  any  point 

in  the  body  of  the  gas  is  determined  by  the  simple  wave  repre- 
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sentation.  In  the  "small  rate"  limit,  <<  ,  the  non¬ 

linear  functional  equation  which  determines  the  signal  in  the 
periodic  state  can  be  reduced  to  a  nonlinear  ordinary  differen¬ 
tial  equation.  Further,  in  this  limit,  the  distortion  of  a 
signal  in  one  travel  time  is  negligible.  The  final  differ2''*ial 
equation  only  describes  the  signal  at  a  boundary  in  the  periodic 


State,  and  cannot  describe  the  cumulative  distortion  which  led 
to  this  state.  It  is  evident  that,  for  a  given  piston  amplitude, 
the  small  rate  limit  restricts  the  allowed  values  of  the  applied 
frequency.  In  fact,  for  some  experimentally  used  values  of  the 
piston  amplitude  the  restriction  is  surprisingly  strong. 

There  are  two  basic  phenomena  in  the  model  used  here; 
shocks  due  to  nonlinearity  and  damping  which  can  prevent  shocks. 
For  a  purely  inviscid  model,  linear  theory  predicts  ar,  unlounded 
amplitude  in  the  periodic  state  for  certain  discrete  (resonant) 
frequencies.  On  the  other  hand,  nonlinear  inviscid  theory  pre¬ 
dicts  a  bounded  signal,  which  contf’^ns  shocks,  in  a  band  about 
the  resonant  frequencies.  Since  shocks  act  as  a  dissipative 
mechanism  they  allow  a  balance  of  energy.  This  role  of  non¬ 
linearity  seems  to  be  well-understood.  To  date,  the  most  com¬ 
prehensive  investigation  of  the  effect  of  damping  on  resonant 
motions  is  due  to  Chester  (1964) .  He  investigated  the  effects 
of  compressive  viscosity  and  boundary  layer  friction  on  the 
motions.  His  conclusion  was  that  the  effect  of  the  former  is, 
in  most  cases,  negligible  except  in  the  interior  of  a  shock.  He 
concurred  with  Betchov  (1958)  that  damping  could  have  a  signifi¬ 
cant  effect,  and  was  led  in  his  analysis  of  boundary  layer 
friction  to  raise  the  question  of  the  existence  of  a  finite 
critical  value  of  damping  which  would  ensure  a  continuous  motion. 
He  concluded  that  it  was  an  open  question.  We  iiiuroduce  three 
different  types  of  damping  which  can  be  treated  within  the  same 
theoretical  framework.  These  are  damping  due  to  radiation  of 
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energy  from  one  end  of  the  tube,  rate  dependence  of  the  gas, 
and  boundary  layer  friction.  We  define  a  lumped  damping  co¬ 
efficient  which  incorporates  the  effects  of  all  three,  and 
show  that  these  "lower  order"  damping  mechanisms  can  prevent 
the  occurrence  of  shocks  in  the  flow,  i.e.,  for  a  given  piston 
motion  there  is  a  critical  level  of  damping  above  which  the 
gas  motion  is  continuous.  B.  Sturtevant,  at  the  California 
Institute  of  Technology,  has  carried  out  experiments  in  which 
a  hole  is  made  in  the  closed  end  of  the  tube.  He  found  that 
for  the  particular  conditions  of  an  experiment  there  is  a  criti¬ 
cal  ratio  of  the  area  of  the  hole  to  the  area  of  the  tube  end 
at  which  shocks  disappear.  This  points  to  the  existence  of  a 
critical  amount  of  damping  which  ensures  a  shockless  motion. 

In  section  3  there  is  a  complete  analysis  of  the  resonant 
band  for  both  the  inviscid  and  damped  cases.  This  is  achieved 
by  an  examination  of  the  integral  curves  of  the  governing  dif¬ 
ferential  equation,  using  a  condition  on  the  mean  of  the  flow 
to  fix  the  shock  position.  Explicit  analytical  results  are 
given  for  the  inviscid  case;  some  qualitative  results  for  the 
damped  motion  are  found  analytically  while  quantitative  results 
are  determined  numerically.  In  section  5  there  is  a  comparison 
between  theory  and  experiment.  The  introduction  of  damping 
improves  the  agreement.  It  suggests  that  if  one  is  interested 
in  such  gross  features  of  the  flow  as  the  maximum  or  minimum 
pressure,  or  the  shock  strength,  then  the  lumped  damping  pro¬ 
vides  an  adequate  description. 
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2.  Formulation 

A  column  o£  gas,  of  length  L  in  some  reference  (equi¬ 
librium)  state,  is  contained  in  a  pipe.  One  end  of  the  pipe 
is  closed  while  at  the  other  end  there  is  an  oscillating  piston. 
If  pressure  and  density  are  measured  from  their  values  in  the 
reference  state  (p^.p^)  with  the  associated  sound  speed  , 
then  in  terms  of  the  nondimensional  variables  (a  u,p  a  p,p  p) 
and  (Lx,  La”^t)  the  governing  equations  in  Lagrangian  form 
are 

[(l+e)-^],^  -  u,^  =  0  ,  (2.1) 

and 

u,^  +  ~  ^  *  C2.2) 

where  e(=p-l)  is  the  condensation,  yp  the  excess  pressure 
ratio  and  u  the  nondimensional  particle  velocity.  The  equation 
of  state  of  a  polytropic  gas  in  these  variables  is 

YP  =  (1+e)'^  -  1  .  (2.3) 

The  end  x=0  is  considered  to  be  'near  rigid*  in  the  sense  that 
we  allow  for  the  possibility  of  radiation  of  energy  through  this 
end  of  the  tube,  but  do  not  consider  the  case  wh^n  it  is  open. 

A  boundary  conditicn  of  this  nr-iture  has  been  discussed  by  Mortell 
and  Varley  (1570) .  Across  the  boundary  “t  x=0  both  pressure 
and  velocity  are  continuous,  and  so  the  disturbance  must  be  com¬ 
patible  with  the  homogeneous  boundary  condition 


w 


uro,t) 


•  ip(0,t)  , 


(2.4) 
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where  i/y  (^0)  is  the  impedance  of  the  interface.  Here  the 
essential  assumption  is  that  the  disturbance  outside  the  tube 
is  generated  by  the  passage  of  a  simple  wave.  Note  that  i-0 
corresponds  to  a  rigid  end  and  to  an  open  end.  We  examine 

the  small  amplitude,  time-periodic  re-sponse  of  the  gas,  governed 
by  equations  (2.1) -(2.3),  to  the  boundary  condition  (2.4)  at 
x=0  and  a  periodic  piston  displacement  at  x=l  of  ;;he  form 
eh(ut}  .  The  amplitude  of  the  displacement  is  e(<<l)  ,  and 
the  period  of  h  is  normalized  so  that  h(y+l)  =  h(y)  .  Then 
the  piston  velocity  at  x=l  is 

uf.  ,'}  *  ewh’(a)t)  =  H(a3t)  .  (2.5) 

Since  h  i.  periodic,  integration  of  (4)  yields 


hquations  (2.1) -(2.3)  are  nonlinear  and  admit  discontinuous 

solutions.  However,  it  has  been  shown  by  Mortell  and  Seymour 

(i972b)  that  for  time-periodic  motions,  be  they  continuous  or 

discontinuous,  the  mean  pressure  and  velocity  do  not  vary  from 

particle  to  particle.  By  choosing  as  the  reference  pressure, 

p^  ,  the  constant  mean  of  the  periodic  state,  conditions  (2,4) 

and  (2.6)  imply  that  the  means  of  u  and  p  are  zero.  The 

actU'».!  value  of  p  can  be  determined  only  from  an  initial 

o 

valur  problem. 

2 . 1  I-qgation  for  the  Periodic  Motion 

A  representation  derived  by  Mortell  and  Varley  (1970)  is 
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used  to  reduce  the  nonlinear  boundary  value  problem  defined  by 
(2.1) -(2.5)  to  a  nonlinear  difference,  or  functional,  equation. 
For  a  more  restricted  class  of  problems,  thi.>  functional  equa¬ 
tion  may  be  further  reduced  to  ?,  nonlinear  ordinary  differential 
equation  which  determines  the  periodic  motion  of  the  gas. 


It  is  convenient  to  reformulate  equations  (2.1) -(2.3)  in 
terms  of  the  Riemann  invariants  and  characteristic  curves  of 
the  system.  Upon  defining 


c(e) 

where 


a(s)(l+s)''^  ds  =  e[l  +  (|M-l)e  •*  0(e^)] 

0 


a^(e)  =  (1+e)^  ^  =  1  +  2Me  +  O(e^)  , 

and  M=^(y+1)  ,  equations  (2.1)-(2.3)  define  the  Riemann  in¬ 
variants 

2f(3)  =  u-c  =  u-p+  O(e^)  (2.7) 

and 

-  2g(a)  =  u+c  =  0+^+  O(e^)  .  (2.8) 


The  associated  characteristics  are  given  by 

*■  a(e)  and  ^  =  -a(e)  .  (2.9) 

a  3 

When  only  one  component  of  the  motion  is  excited,  equations 

(2.7) -(2.9)  admit  two  exact  solutions,  simple  waves,  which 

correspond  to  f=constant  and  g=con?tant.  When  both  components 

of  the  motion  are  excited  there  is  in  general  an  interaction 


dx 

3t 


-8- 


between  a-waves,  moving  to  the  right,  and  8-waves,  moving  to 
the  left.  However,  it  has  been  shown  by  Mortell  and  Varley 
(1970)  that  to  first  order,  in  the  limit  of  small  amplitudes, 
the  waves  do  not  interact  as  they  pass  through  each  other  in 
the  body  of  the  gas.  By  this  is  meant  that  to  first  order  the 
trajectory  of  an  a-wave  is  determined  only  by  the  signal  it 
carries  and  is  not  influenced  by  the  8-w'  /es  through  which  it 
passes.  Thus  the  motion  of  the  gas  may  be  represented  as  the 
s;  t 'srposition  of  two  noninteracting  simple  waves.  Then  equations 
(2.7)  and  (2.8)  imply  that  to  first  order 

e  =  p  =  -  f(8)  -  g(a)  and  u  =  f(B)  -  g(a)  ,  (2.10) 

while  (2.9)  integrate  to  give 

^  =  t  -  X  -  Mxg(a)  and  ^“t+x-l+  M(x-l)f(8)  ,  (2.11) 

where  we  have  parametriiied  a  and  8  by  a=a)t  on  x“0  and 

8®a)t  on  x=l  .  Upon  using  the  boundary  conditions  (2.4)  and 

(2.5),  g  is  eliminated  from  (2.10)  and  (2.11)  to  yield  the 
nonlinear  functional  difference  equation  to  determine  tiu-  signal, 
f  ,  on  the  boundary  x=l  : 

f(n)  -  kf(s)  =  H(ri)  (2.12) 

where 

n  =  s  +  2w  +  a)M(l+k)f(s)  .  (2.13) 

In  (2.12)  and  (2.13),  k  =  is  the  reflection  coefficient 

at  x=0  ,  where  g  is  related  to  f  by 
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g(0)  =«  k£((j))  ,  (2.14) 

and 

0  =  (J)  +  0)  +  wMfC({))  .  (2.15) 

The  governing  differential  equations  (2.1) -(2.3)  and  the 
boundary  conditions  (2.4)  and  (2.5)  have  been  reduced  in  the 
small  amplitude  limit  to  the  functional  difference  equation 
(2.12)  and  (2.13).  We  nov/  seek  solutions  to  (2.12)  and  (2.13) 
which,  like  the  piston  motion,  have  unit  period.  Further,  as 
a  consequence  of  the  representations  (2.10)  and  (2.11),  the 
boundary  conditions  (2.4)  and  (2.5)  and  the  fact  that  u  and  p 
have  zero  mean  over  any  period,  f  and  g  must  satisfy 

[^f(s)ds  =  [^g(s)ds  =  0  .  (2.16) 

^0  0 

Since  M  is  the  ratio  of  second  order  to  first  order  elastic 
constants,  linear  theory  is  recovered  from  (2.12)  and  (2.13) 
by  formally  setting  M=0  to  yield 

f(n)  -  kf(n-2w)  =  H(n)  .  (2.17) 

I'fhen  k=l  there  are  no  solutions  of  (2.17)  with  unit  period; 
that  is ,  when 

,  a)  =  w  =in,  n  =  l,2,3,...  .  (2.18) 

n  d* 

These  are  the  linear  resonant  frequencies.  Ultimately  we  will 
consider  the  time-periodic  response  of  the  system  to  tj?quencies 
near  to  those  defined  by  (2.18)  and  consequently  define 

u  »  u  (l+S)  '  ( I  6  I <^)  • 

n  ^  ‘  '  n 


(2.19) 
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Then  in  terms  of 

FCy)  =  f(y)  +  ^  .  (2.20) 

where 

b  =  ^(l+k)Cl+6)  (=0(1))  ,  (2.21) 

(2.12)  and  (2.13)  become 


F(n)  -  kFCs)  =  GCn)  , 


(2.22) 


and 


n  =  s  +  n  +  nbF(s)  ,  (2.23) 

where 

G(n)  =  ^  +  H(n)  and  m  =  1  -  k  .  (2.24) 

Definition  (2.20)  now  implies  that  the  zero  mean  condition, 
(2.16)  on  f  is  replaced  b/ 

f%(s)ds  =  i  .  (2,25) 

•'  0 

The  approximations  used  to  derive  equations  (2 . 22) - (2. 24) 
are  the  small  amplitude  assumption,  |f|  ,  |g|  <<  1  ,  and  the 
fact  that  the  impedance  of  the  interface  at  x=0  is  near  zero, 
sc  that  0  k  £  1  .  This  latter  assumption  is  required  since 
whenever  1  +  k  =  0(e)  the  next  correction  to  the  character- 
istics  (at  order  e  )  is  no  longer  negligible.  For  example, 
if  the  end  of  the  tube  is  open  (k=-l)  tl^  nonlinear  approxima¬ 
tion  (2.11)  to  the  characteristics  leads  to  a  linear  difference 
equation  which  has  no  bounded  periodic  solution  at  a  resonant 
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frequency.  Consequently,  for  a  problem  involving  an  open  end 
(or  "nearly*'  open)  the  approximations,  (2.11),  for  a  and  g 
must  be  improved.  In  fact,  the  approximations  for  a  and  3 
must  contain  terms  in  f  and  g  so  that  now  the  motion  in 
the  tube  is  determined  by  the  cubic  term  in  the  equation  of 
state,  with  a  resulting  amplitude  at  0(e'^^)  .  Note  further 
that  the  difference  equation  (2 . 22) - (2 . 24)  together  with  (2.20) 
determine  the  shape  of  the  signal  function,  £  ,  only  on  x=l  . 
The  velocity  and  pressure,  u  and  p  ,  are  subsequently  calcu¬ 
lated  at  any  particle  x  in  the  tube  from  the  representations 
(2.10)  and  (2.11).  This  is  particularly  important  when  there 
is  significant  distortion  of  a  wave  in  one  cravel  time. 

If  we  now  make  the  additional  small  rate  assumption 

|nF*|  «  1  , 


equation  (2.23)  implies 

F(s)  =  F(n-n)  -  nbF(n-n)F'(n-n) [l+0(nF')]  .  (2.26) 

Upon  using  (2.26),  and  since  we  seek  solutions  with  unit  period, 
the  difference  equation  (2.22) -(2.24)  can  then  be  approximated 
by  the  nonlinear  ordinary  differential  equation 

vF(n)F’(n)  +  pF(n)  =  G(n)  ,  0<ri<l  .  (2.27) 

together  with 

F(n+1)  =  F(n)  ,  (2.28) 


where  v  =  nbk  .  The  small  rate  condition  and  the  definition  of 
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F  give  the  further  restriction 

|6|  «  1/n  .  (2.29) 

Thus  equation  (2.27)  is  valid  in  the  small  rate  limit  only  for 
periodic  motions  at  frequencies  in  the  neighborhood  of  linear 
resonant  frequencies.  In  contrast  to  this,  the  difference 
equation  (2.22) -(2.24)  was  derived  with  no  restriction  on  the 
applied  rates;  it  is  valid  for  non-periodic  phenomena  (see 
Mortell  and  Varley  (1970),  Mortell  and  Seymour  (1972a, b)),  and 
the  applied  frequency  is  not  restricted  to  lie  near  a  linear 
resonant  frequency. 

Since  the  small  amplitude  restriction  requires  that  (f(<<l  , 
the  small  rate  condition  implies  that  the  differential  equation 
(2.27)  is  a  good  approximation  to  the  difference  equation  only 
when  n  is  at  most  0(1)  as  lf|-*-0  .  In  addition,  the  repre¬ 
sentations  (2.11)  for  the  characteristics  imply  that  in  this 
limit  there  will  be  no  appreciable  distortion  of  the  waveform 
in  one  period,  since 

6  =  w(t+x-l)  +  0(nf)  . 

(Of  course  there  will  be  a  cumulative  distortion  of  the  signal 
until  the  periodic  motion  has  been  set  up.  This  is  not  described 
by  (2.27)).  Hence,  in  the  small  rate  limit,  (2.10)  and  (2.11) 
can  be  replaced  by  the  linear  acoustic  representation 


-p  »  f(w[t+x-l])  +  g(w[t-x])  and  u  «  f(u[t+x-l])  -  g(u[t-x])  , 

(2.30) 
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so  that  when  the  periodic  motion  has  evolved,  nonlinearity  is 
of  primary  importance  in  determining  the  shape  of  the  signal 
functions,  but  it  is  of  secondary  importance  in  determining 
how  the  signals  propagate. 

We  point  out  that  the  small  rate  condition  |nF*!  <<  1  , 
which  is  necessary  for  the  validity  of  both  the  differential 
equation  (2.27)  and  the  representation  (2.30),  is  quite  re¬ 
strictive  for  usual  experimental  values  of  the  parameters.  The 
results  of  section  3  show  that  for  e  =  .0147  ,  used  by  Stur- 
tevant  (1972),  this  restriction  implies  n<<2  ;  for  t  =  .0018  , 
used  by  Saenger  and  Hudson  (1960) ,  the  restriction  yields  n<<6  . 
In  a  sequel  to  this  paper  we  shall  examine  periodic  motions  with 
no  restrictions  on  the  rates  when  some  of  the  ideas  introduced 
here  are  used  to  analyze  the  functional  difference  equation, 
(2.12)  and  (2.13),  directly. 

It  may  be  of  interest  to  note  that  equation  (2.27)  arises 
in  other  physical  situations.  It  is  a  generalization  of  the 
equation  which  describes  the  motion  of  a  viscously  damped  pendu¬ 
lum  under  a  constant  external  moment,  and  also  occurs  in  the 
study  of  the  pull-out  torque  of  a  synchronous  motor,  see  Stoker 
(1950)  or  Minorsky  (1962).  In  Appendix  I  we  show  how  equation 
(2.27)  may  be  derived  from  the  governing  equations  (2.1) -(2.6) 
by  a  regular  perturbation  procedure,  and  point  out  an  extra 
limitation  on  p  imposed  by  that  procedure.  The  analysis  of 
equation  (2.27),  under  the  restriction  (2.25),  and  the  physical 
interpretation  of  the  results  constitute  the  remainder  of  this 


paper. 
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3.  Determination  of  the  Periodic  Signal 

Here  we  analyze  the  integral  curves  of  the  differential 
equation  (2.27)  for  various  ranges  of  the  parameters  y  and  6 
and  then  use  these  to  construct  the  signal,  F  ,  of  the  per¬ 
iodic  motion.  When  we  have  constructed  F  over  one  period, 
possibly  by  a  composition  of  integral  curves,  it  is  continued 
periodically  by  (2.28).  The  signal  function,  F  ,  may  then 
have  discontinuities  representing  a  time  periodic  motion  in  the 
pipe  containing  shocks.  The  discontinuities  in  F  arise  in 
satisfying  the  mean  condition  (2.25).  Acoustic  theory  allows 
discontinuities  of  either  compression  or  rarefaction  with  no 
restriction  on  .their  strengths.  However,  to  be  physically 
acceptable  a  jump  in  a  gas  must  be  compressive.  Here  we  con¬ 
sider  only  piston  velocities  which  have  three  zeroes  over  one 
period,  and  then  discontinuous  solutions  of  (2.27),  which  sat¬ 
isfy  the  mean  condition  (2.25),  additionally  satisfy  the  appro¬ 
priate  weak  shock  relations.  This  is  not  strictly  necessary 
within  the  acoustic  approximation. 

If  S(x)  is  the  arrival  time  at  x  of  a  weak  shock 
traveling  in  the  negative  x  direction  and  0~(x)  are  the 
wavelets  immediately  ahead  of  and  behind  the  shock,  then  the 
weak  shock  relations  imply  that 

^  -  1  -  ♦  fC6-)]  .  (3.1) 

A  similar  relation  gives  the  speed  of  shocks  moving  to  the 
right.  However,  since  in  the  periodic  state  there  is  negligible 
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distortion,  0"  are  independent  of  x  and  (3.1)  can  be  in¬ 
tegrated  to  give  the  travel  time  of  a  shock  from  x=l  to  x=0  . 
The  boundary  condition  (2,14)  then  implies  that  the  total  travel 
time  for  the  shock  to  return  to  x=i  is 

T  =  2  +  |M(l+k)[f(0'")  +  f(B-)]  .  (3.2) 

Since  the  periodicity  requirement  is  that  ^  ^  ,  (3.2)  implies 

that 

f(B-")  +  f(B-)  =  -  ^  •  (3.3) 

On  using  (2.20),  the  definition  of  F  ,  (3.3)  becomes 

F(B'^)  +  F(B-)  =  0  .  (3.4) 

The  condition  that  only  compressive  shocks  are  allowed  then 
requires  that  F(B’’')  >  0  . 

On  the  other  hand,  integration  of  (2.27)  over  one  period, 
assuming  a  discontinuity  at  n^B  yields 

|v[f2(B'')  -  f2(B-)]  +  u  [^F(s)ds  =  [^G(s)ds  .  (3.5) 

Conditions  (2.6),  (2.24)  and  (2.25)  then  imply  that  (3.4)  must 
hold  at  the  discontinuity.  Thus  a  solution  of  (2.27),  contain¬ 
ing  a  discontinuity,  which  satisfies  the  mean  condition  (2.25) 
and  the  restriction  F(B'*’)  >  0  v^/ill  necessarily  satisfy  the 
weak  shock  relations.  Thus  a  shock  is  fitted  into  the  solution 
by  satisfying  the  mean  condition.  This  analysis  is  for  one 
shock  per  period  of  the  piston,  which  gives  n  shocks  in  the 
tube  at  any  time. 
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3 cl  Special,  case  of  ao  damping^  y“0 

When  the  boundary  at  x“0  is  rigid  and  there  is  no  radia¬ 
tion  o£  energy  through  it,  i^y^O  and  equation  (2.27)  is  greatly 
simplified.  It  can  then  be  integrated  cofljpletely  and  the  signal 
function,  the  width  of  the  resonant  band  and  the  shock  strength 
for  these  frequencies  can  be  determined  arf'’,lytically  in  terms 
of  the  parameters  of  the  problem.  Further,  the  transition  from 
a  discontinuous  motion  inside  the  band  to  a  continuous  one  out¬ 
side  is  exhibited  explicitly. 

We  wish  to  distinguish  between  the  integral  curves  of  equa¬ 
tion  (2.27)  and  the  signal  F  which  must  additionally  satisfy 
the  mean  condition  (2.25).  An  integral  curve  is  denoted  by  Z(ti)  . 
Notice  that  while  F  is  defined  only  for  05n£l  and  is  then  con¬ 
tinued  periodically,  if  F  is  continuous  it  must  coincide  with  an 
integral  curve  Z  which  is  both  continuous  and  periodic  for  -oo<n<oo 
Conversely,  a  continuous,  periodic  integral  curve  Z  with  unit 
period  which  satisfies  the  mean  condition  (2.25)  is  the  required 
signal  function  F  .  When  such  an  integral  curve  exists  it  is 
unique.  When,  for  a  particular  frequency,  no  such  curve  exists, 

F  is  discontinuous  and  is  composed  of  distinct  integral  curves. 

For  the  case  ij=0  ,  the  appropriate  differential  equation  is 

nbZ(n)Z’(n)  =  H(n)  ,  (3.6) 

There  is  no  loss  of  generality  in  choosing  the  origin  so  that 
H(0)  =  H(n^)  =  H(l)  =  0  ,  where  0<nj^<l  ,  with  H’(0)  =  H'(l)  >  0 
and  H’(n^)  <  0  .  Then  in  the  (n,Z)  plane  the  points 
\  =  (0,0)  ,  Ag  =  (1,0)  and  B  ^  are  isolated  singular 
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pointf^.  and  Ag  are  saddle  points,  while  is  a  center. 

The  two  separatices  through  A^  are  given  ly 

Z-Cn)  =  -[^  |'’H(s)ds]^  ,  (3.7) 

while  the  zero  mean  condition  (5)  implies  that  Z"(l)  =0  . 

Thus  the  separatices  connect  the  two  saddle  points  A^  ,  A^  . 

The  other  integral  curves  which  are  defined  for  all  n  are 
given  by 

zCn)  =  |\cs)ds  (3.8) 

where  C=Z^(0)  .  These  solutions  are  periodic,  with  imit 
period,  for  all  n  • 


Since  Z'^(r))  >  0  for  0<n<l  ,  any  solution  Z(n)  with 
Z(0)  >  0  is  periodic  in  n  with  Z(ri)  >  Z’‘'(n)  and  therefore 
Z  satisfies 


Z(s)ds  >  Z'^(s)ds  . 


(3.9) 


Consequently  the  mean  condition  (2.25)  implies  that  for  an 


applied  frequency  w  =  ^(1+6)  such  that 


i  > 


rl 


Z'^(s)d- 


(3.10) 


there  exists  a  unique,  continuous ,  periodic  solution 


z.Cn) 


5F 


^’^H(s)ds  +  C(6)]*^  =  A(n)  .  (3.11) 


The  positive  constant  C(6)  is  choosen  so  that  Zg(n)  satisfies 
the  mean  condition  (2.25).  A  similar  analysis  shows  there  are 


-18- 


also  continuous  periodic  solutions, 

Zg(n)  =  -  A(n)  for  ^  <  f  Z-(s}ds  <  0  .  (3.12) 

Thus  for  the  range  of  frequencies  defined  by  (3.10)  and  (3.12) 
the  signal  function  F  is  given  by 


F(n)  =  Z^(n) 


(3.13) 


and  is  continuous  and  periodic. 


For  frequencies  such  that 


[^Z-(s)ds  <  ^  <  [^Z+(s)ds  (3.14) 

no  single  integral  curve  wrill  satisfy  the  mean  condition  (2.25) 
and  the  signal  function  will  necessarily  be  discontinuous .  The 
shock  condition  together  with  the  fact  that  only  comp.-essive 
shocks  are  allowed  then  implies  that  the  signal  function  F  can 
only  be  constructed  from  the  separatices  Z“(ri)  with  Mist  one 
shock  per  period.  The  position  of  the  shock  at  n=n„  is  chosen 
to  satisfy  the  mean  condition.  The  signal  function  F  is  then 
given  by 


and 


F(n) 


ns<nii 


F(n+1)  =■  F(n)  , 


(5.15) 


where  condition  (2,5)  implies  that 


|{h(n)-h(0)} 


(3J.6) 
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The  range  of  frequencies,  defined  by  (3.14),  for  which  the 
signal  is  discontinuous,  is  called  the  resonant  band.  If  (3.14) 
is  solved  for  6  the  resonant  band  is  given  e.Aplicitly  by 

6-  <  6  <  6+ 


where 


± (eM)^h 
l+(eM)^h 


(3.17) 


and  h  =  [h(s)-h(0)]^  ds  ,  Notice  that  the  amplitude  of  the 

response  of  the  gas  to  an  applied  sig^ial  of  0(e)  is  O(e^) 
and  that  the  width  of  the  issonant  band 


l-eMii^ 


is  aJ.'.o  O(c^)  . 


The  above  results  are  particularly  simple  for  the  important 
special  forcing  function 


h(n)  =  -  cos2Trn  . 


Then 


Z-(r.) 


*  (r) 


simrn 


(3.18) 


(3.19) 


and  ii  =  Z/I/tt  .  In  addition  the  shock  strengch  can  be  found 
explicitly  in  terms  of  6  ; 


Shock  strength  •*  Z‘^(n  )  -  Z"(n  ) 

s  s 

=  2(^)^(1-4.2(6))*^  , 


(3.20) 


where 


2(l  +  6)(2eM)' 


and  cos  (frn  ) 
s 


(3.21) 
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Thus  for  a  given  6  within  the  resonant  band,  the  shock  strength 
is  given  by  (3.20)  and  the  position  of  the  shock  by  the  second 
of  (3.21).  On  using  (3.17),  it  follows  that 

4,(6+)  =  -  4,(6-)  =  -  1  , 

and  so,  by  (3.20),  the  shock  strength  tends  to  zero  as  6-»-6“  . 
Further,  as  6-+-6‘'’  ,  Hg-^l  ,  while  as  6->6“  ,  Hg-^-O  .  The  limiting 

•4*  4’ 

solutions,  when  6=6“  ,  are  given  by  F(n)  =  3"(n)  ,  0<n£l  . 

Then  the  signal  F  is  continuous  but  has  a  discontinuous  slope 
at  ri=0  and  n=l  •  The  resonant  band  is  not  symmetrically 
situated  about  the  linear  resonant  frequencies  o)=co  since 
|<S^|  >  |6~1  ,  by  (3.17).  Note,  further,  that  F'(n^)  =  0  and 
thus  for  the  case  h  "  -cos2TTn  and  w=a)  ,  n  =  =  4  »  i-e., 

the  maximum  pressure  equals  the  pressure  immediately  ahead  of 
the  shock,  and  the  pressure  immediately  behind  the  shock  is  the 
minimum  pressure. 

The  main  point  of  the  analysis  given  above  Cor  the  undamped 
case  is  that  all  the  results  are  obtained  from  the  simple  equa¬ 
tion  (3.6)  together  with  the  mean  condition  (2.25).  By  this  we 
mean  that  for  any  frequency  w  such  that  <<  1  (either 

inside  or  outside  the  resonant  band)  the  signal  function,  F  , 
is  determined  by  (3.6)  and  (2.25).  This  is  a  consequence  of  the 
translation  (2.20).  The  advantage  of  this  approach  is  that  it 
does  not  rely  on  having  an  integral  of  the  governing  equations 
to  construct  the  resonant  band,  and  hence  can  be  generalized  to 
the  dissipative  case,  yfO  .  The  resonant  band  for  the  case 
y=0  ,  h  =  -cos2Trn  has  been  treated  previously  by  Chester  (1964). 
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His  analysis  of  the  band  involves  finding,  by  an  ad  hoc  pro¬ 
cedure,  an  equation  which  is  uniformly  valid  in  the  frequency 
parameter. 

3.2  <  y  <  1 

When  y  is  nonzero  the  positions  of  the  singular  points 

in  the  (HiZ)  plane  now  depend  on  both  y  and  6  and  is 

no  longer  a  center.  We  consider  the  variations  in  (y,6)  in 

two  parts.  Here  we  show  that  for  a  given  forcing  function  ll(n) 

there  is  a  critical  amount  of  damping,  y=y^  ,  such  that  for 

y>y^  the  signal  function  is  continue  us  for  all  frequencies. 

In  section  3.3  we  fix  y<y  and  consider  variations  in  which 

c 

will  define  the  resonant  band. 

If  we  assume  that  the  zeroes  of  H(ri)  satisfy  tne  conditions 
described  in  3.1,  the  singular  points  of  the  equation 

vZ(n)Z'(n)  +  yZ(n)  =  GCn.)  ,  (3.22) 

where 

G(n)  =  H(n)  +  P  ^  .  (3.23) 

are  the  points  (9^,0)  ,  i=0,l,2,  such  that  G(6^)=0  .  (The 
periodicity  of  H  ensures  that  6,  =  9^  +  1  .)  We  label  them 
A^(y,6)  ,  A2(y,6)  and  Bj^(y,6)  where  A^(y,0)  =  (0,0)  , 

A2(y>0)  =  (1,0)  and  B^(y,0)  =  (t^^,C)  .  Labeling  them  in  this 
way  is  consistent  with  the  notation  of  the  previous  sectit r  and 
ensures  that  A  and  A^  avs  again  saddle  points.  The  sty  ir- 
atices  through  A^  (which  wj  denote  by  ZT(n)  )  have  slopes 
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A-C0^}  =  |-y±[y2+4vH'(0.)]^y2v  ,  (3.24) 

where  A‘*'C0^)  >  0  >  X“(6^)  ,  since  H'C6^)  >  0  ,  for  i=0,2  . 

The  slopes  at  are  also  given  by  (3.24),  where  H*(0^)  <  0  . 

When  >  0  >  where 

I(n)  =  +  4vH'(n)  ,  (3.25) 

is  a  node,  while  if  1(0^)  <0  B^  is  a  focal  point.  Thus 
for  a  given  forcing  function  H  (which  is  0(c)  )  Bjl  will  be 
a  node  if  there  is  sufficient  damping  in  the  system.  It  may  then 
be  possible  to  construct  a  continuous  solution  passing  through 
,  B^  and  for  any  value  of  6  .  Obviously  the  nodal 

condition  1(0^^)  >  0  is  necessary  for  the  existence  of  such  a 
solution,  however  it  may  not  be  sufficient.  Since  the  distortion 
of  the  signal,  and  possible  shock  formation,  depends  on  the 
amplitude  of  H'(n)  ,  one  can  expect  the  condition  ensuring  the 
existenc-:'  of  a  continuous  solution  to  depend  on  a  global  property 
of  H ' (n)  .  In  fact 

y^  >  y^  =  max  [-4vH'(n)]  >  0  (3.26) 

n 

is  a  sufficient  condition  for  the  existence  of  a  continuous, 
periodic  solution  at  all  frequencies.  The  proof  of  this  result 
is  given  in  Appendix  2.  Thus  for  y>y^  , 

F(n)  =  Z^(n)  ,  0<n<i  ,  (3.27) 

which  is  continuous  and,  by  ''3.5),  satisfies  the  mean  condition 
(2.2,5),  (see  figure  [4])* 
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3.3  0  <  M  <  y 

c 

In  this  case  there  is  not  enough  damping  to  produce  a 
shockless  solution  for  all  frequencies.  Thus  for  6=0 

Z^(0^)  >  0  ,  Z“(03^)  <  0 

while 

where  O<n2^0j^<0Q‘^l  (see  figure  (  1  ])  .  Hence  the  separatices 

do  not  connect  the  saddle  points  A  and  Ap  .  As  6  is  in- 

creased  through  the  resonant  band,  so  that  G(s)ds  >  0  ,  there 

'  o 

exists  a  unique  frequency,  given  by  6=6'*'  ,  such  that 

z+(n)  S  z+(n)  ,  Qpinieg  . 

That  is,  for  6=6'*'  ,  the  positive  separatrix  connects  the  saddle 

points  A^  and  ,  (see  figure  [2  ]).  Further,  for  6>5'^ 

there  exists  a  unique,  continuous  periodic  solution  Z=Zg(n)>0  , 

(see  figure  [3]).  Similarly  there  exists  6=6“<0  for  which 

Z~=Z“  and  such  that  when  6<6“  there  is  a  unique ,  continuous 

periodic  solution  Z=Zg(n)<0  .  These  results  can  be  inferred 

from  the  results  of  Amerio  (1949,1950).  Whereas  for  the  case 

p=0  explicit  values  have  been  given  for  6"  (see  equation 

(3.17)),  when  0<m<m^  this  is  not  possible.  However,  for  a 

particular  forcing  function,  6"  Can  easily  be  found  numerically 

by  varying  6  until  a  solution  is  found  such  that  Z^(02)=O  or 

Z"(0  )=0  .  Since  these  limiting  solutions  are  continuous  they 
2  o 

satisfy  the  mean  condition  (2.25)  and  hence  we  can  give  the  im- 
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plicit  conditions  for  the  edge  of  the  band  as 

'  mF)-  ^  . 

0  O 

There  have  been  several  attempts  to  obtain  analytical  bounds  on 
6~  .  Hayes  (1953)  and  Bdhm  (1953)  found  bounds  for  h=-cos2Trn 
while  Lillo  and  Seifert  (1955)  used  similar  techniques  to  find 
bounds  for  a  general  forcing  function.  Further  reference  can  be 
found  in  Sansone  and  Conti  (1964) 

By  equation  (3.5)  the  unique  continuous  solutions  Z^(ri) 
automatically  satisfy  the  mean  condition  (2.25)  and  hence  for 
or  6<6“ 

F(n)  =  z^(n)  ,  •  (3.28) 

V/e  contrast  the  case  of  no  damping,  y=0  ,  when  the  continuous 
periodic  solutions  Z^(n)  were  not  unique,  but  the  mean  condition 
uniquely  determined  the  signal  function  F  . 

When  6~<6<6’^  there  are  no  continuous  periodic  solutions 
of  (3.22).  Again,  like  the  case  y=0  ,  we  construct  the  signal 
function  F  by  a  composition  of  integral  curves,  (see  figure 
[l]).  The  discontinuous  signal  function  F  must  satisfy  both 
the  mean  condition  (2.25)  and  the  weak  shock  condition  (3.4). 
However  it  has  been  shown  that  if  the  mean  condition  is  satisfied 
the  shock  condition  is  automatically  satisfied.  The  condition 
that  a  shock  is  compressive  then  implies  that  we  choose 
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F(n) 


ng<nie2 


(3.29) 


where  is  the  position  o£  the  shock.  It  is  shown  in 

Appendix  2  that  it  is  always  possible  to  choose  an  n_  to  com- 
bine  Z^(n)  and  Z“(n)  so  that  the  mean  condition  (2.25)  is 
satisfied.  Hence  F  ,  as  given  by  (3.29),  is  the  required 
signal  function.  By  equation  (3.22),  F'=0  at 

where  F=G/p  .  Then,  in  general,  when  the  maximum  pres¬ 

sure  exceeds  the  pressure  immediately  ahead  of  the  shock  and 
the  pressure  immediately  behind  the  shock  exceeds  the  minimum 
pressure. 


It  is  clear  from  the  structure  of  the  integral  curves  for 
the  inviscid  case  that  if  the  piston  frequency  is  an  even  multiple 
of  the  fundamental,  then  a  possible  continuous  solution  is 


with 


F(n)  =  ,  e^<n<02 

F(n+1)  =  -FTnl  , 


which  is  a  "subharmonic"  solution.  From  the  preceeding  analysis, 
this  solution  is  unstable  to  perturbations  in  both  damping  and 
frequency. 

Finally  we  note  that  when  k=0  (y=l)  the  impedances  at  x=0 
are  ’matched’.  Then  there  is  no  reflected  wave,  so  that  by 
(2,14)  giO  .  The  differential  equation  (2.27),  together  with 
(2.20)  and  (2.24)  then  yields  that  on  x=l  , 


f(n)  =  H(n)  for  all  n  . 
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4 .  Lumped  Damping;  Critical  Acceleration  Level 

In  the  preceding  we  have  analyzed  resonant  oscillations 
when  the  only  damping  in  the  system  is  due  to  radiation  of 
energy  auay  from  one  end  of  the  pipe.  Here  we  discuss  two  other 
forms  of  damping  which  fit  into  the  same  theoretical  framework. 
These  are  damping  due  to  internal  dissipation  and  to  wall  fric¬ 
tion.  It  is  shown  how  to  define  a  lumped  damping  coefficient, 
k  ,  which  allows  one  to  incorporate  the  effect  of  the  three 
damping  mechanisms  mentioned  into  an  equation  of  the  form  (2.27). 
The  coefficient  is  referred  to  as  •'lumped"  since  on  the  basis 
of  experiments  it  can  be  used  as  a  measure  of  the  effective 
damping  in  the  system  without  being  able  to  assess  the  effect 
of  the  individual  contributions. 

It  has  been  shown  by  Mortell  and  Seymour  (1972a)  that  the 
representation  (2.10)  and  (2.11)  can  be  extended,  in  the  high 
frequency  limit,  to  include  the  effect  of  internal  dissipation 
of  the  transmitting  media  (specifically,  there,  for  a  viscoelastic 
rod) .  For  a  gas  such  dissipation  would  result  from  the  excita¬ 
tion  of  any  of  the  internal  degrees  of  freedom,  e.g.  vibrational 
excitation  or  molecular  dissociation.  If  it  is  assumed  that  only 
one  rate  dependent  process  is  of  significance  and  that  this  can 
be  represented  by  the  relaxation  variable  a(X,t)  ,  then  the 
rate  of  adjustment  of  a  may  be  described  by 

If  =  '|/(p,e,a)  .  (4.1) 

Ip  will  then  define  a  rate  parameter  or  relaxation  time  ^(>0) 
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proportional  to  ’i'»(jCPo  »®o  ^  periodic  disturbance  is 

considered  of  high  frequency,  or  is  near-frozen,  if  its  period 
is  small  compared  with  x  ,  so  that  LCa^to)”^  <<  x  (note  that 
in  our  variables  LCa^w)"^  represents  the  dimensional  piston 
frequency).  For  such  disturbances  'a^'  should  now  be  inter¬ 
preted  as  the  frozen  rather  than  equilibrium  sound  speed  (more 
details  of  nonlinear  wave  propagation  in  a  relaxing  gas  are 
given  by  Blythe  (1969)).  In  this  limit  a  disturbance  in  the 
gas  can  be  represented  as  two  noninteracting,  modulated  simple 
waves  traveling  in  opposite  directions  (see  Mortell  and  Seymour 
(1972a))  (modulated  simple  waves  in  rate -dependent  media  are 
discussed  in  detail  by  Seymour  and  Varley  (1970)).  The  appro¬ 
priate  representation  corresponding  to  equations  (2.10)  and 
(2.11)  is  then 


e  =  p  =  -f(3)  -g(a)  e-^^  u  =  f(6)  -  g(a)  e""^  , 

(4.2) 

^  =  t  -  X  -  Mg(<x)  (4.3) 

and 

i  =  t  +  X  -  1  +  Mf(3)  d~^(e^^'‘'^^-  1)  ,  (4.4) 

where 


d  =  L/a^x  <<  w  . 


(4.5) 


On  eliminating  g  from  (4.2)  and  (4.3)  through  the  boundary 
conditions  (2.4)  and  (2.5),  the  functional  difference  equation 
(2.12)  and  (2.13)  is  recovered  with  the  parameters  k  and  b 
replaced  by  k  and  b  where 
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ic  =  and  b  =  ^(1+0  (l+ke'*^)  (l-e"^)d"^  .  (4.6) 

In  the  small  rate  limit,  when  a)=0(l)  ,  the  condition  (4.5) 
implies  d<<l  .  Under  these  circumstances  the  procedure  of 
section  2  then  leads  to  the  nonlinear  ordinary  differential 
equation  (2.27)  with  the  parameters  p  and  v  replaced  by 
{i  and  V  where 

ii  =  1  -  k  and  v  =  nSk  .  (4.7) 


The  rate  independent  case  is  recovered  in  the  limit  when 

d->0  .  The  parameter  k  in  (4.6)  consists  of  two  factors;  the 
first,  k  ,  is  the  attenuation  of  the  signal  at  the  interface 
x=0  due  to  transmission  of  energy  into  the  neighboring  medium; 
the  second,  e~’  ,  is  the  attenuation  of  the  signal  over  one 

cycle  due  to  internal  damping.  The  latter  acts  continuously 
throughout  the  gas,  whereas  the  former  only  acts  at  the  inter¬ 
face.  The  role  of  both  in  preventing  shock  formation  is  the 
same . 

In  the  high  frequency  limit,  internal  dissipation  is  a  lower 
order  damping  in  the  sense  that  the  appropriate  linearized  equa¬ 
tion  satisfied  by  u  is 

u.  .  -  u  +  2d  u.  =  0  .  (4.8) 

tt  XX  t 

In  contrast,  a  higher  order  damping,  introduced  for  instance  by 
compressive  viscosity,  in  represented  in  the  linear  case  by 


u 


tt 


XX  o  xxt 
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Higher  order  damping  will  only  structure  a  shock  in  a  resonant 
oscillation  (see  Chester  (1964)),  while  lower  order  damping  can 
prevent  shock  formation  even  for  d<<l  .  A  form  of  damping  which 
is  always  present  to  some  degree  in  a  tube  is  that  due  to  the 
viscous  effects  of  the  boundary  layer  at  the  tube  wall.  Its 
effect  can  be  modeled  as  in  Chester  (1964)  ,  as  a  weighted  in¬ 
tegral  of  F  in  the  equivalent  of  (2.27),  or,  equivalently  by 
adding  a  body  force  term  in  the  momentum  equation.  The  latter 
procedure  was  used  by  Rayleigh  (1945)  for  periodic  oscillations 
in  a  circular  pipe  and  yields  a  lower  order  damping  of  the  form 
given  by  equation  (4.8)  with  the  parameter  'd'  replaced  by 


I  kwL 

la  r^j 


where  k  is  the  kinematic  viscosity  of  the  gas  and  r  is  the 
radius  of  the  tube.  The  effect  is  to  produce  a  damping  coeffi¬ 


cient 


Thus  the  lumped  damping  coefficient  to  account  for  these  three 
damping  mechanisms  takes  the  form  .  The  damping 

coefficient,  k  ,  in  (2.22),  and  the  corresponding  y  in  (2.27) 
can  alv/ays  be  interpreted  in  this  broader  sense.  The  effect  of 
k<l  is  to  decrease  the  shock  strength  in  a  resonant  oscillation. 
Then,  the  effective  damping  in  the  system  can  be  measured  by 
adjusting  k  to  ensure  agreement  between  the  theoretical  and 
experimental  shock  strengths.  A  purely  inviscid  model  predicts 
the  shock  strength  with  an  error  of  about  30%  at  the  fundamental 
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1 inear  resonant  frequency. 

Since  A^Ct)  =  wH*(ut)  represents  the  acceleration  of  the 
piston,  we  can  interpret  condition  (2.26)  as  follows:  for  a 
given  amount  of  damping,  k  ,  there  is  a  critical  acceleration 
level  of  the  piston  such  that  for  applied  accelerations  below 
this  the  motion  of  the  gas  is  shockless  at  all  frequencies.  By 
(3.26),  for  a  known  y  the  gas  motion  is  shockless  provided 
the  piston  acceleration,  satisfies 

|A  (t)|  <  (4.9) 

^  4v  4Mk(l+k) 

on  using  (2.18),  (2.19)  and  (2.21).  In  the  limit  as  the  damping 
tends  to  zero,  k->-l  and  y->-0  which  implies  |Ap(t)|->0  for  a 
shockless  motion.  Thus  when  there  is  no  damping  present  chere 
is  always  a  shock  at  resonance. 

When  i<<l  ,  a<<l  ,  then  (4.9)  reduces  to 

!MAp(t) 1  <  |(i+a)^  (4.10) 

as  the  condition  for  a  shockless  solution  for  a  resonant  forced 
motion.  In  contrast,  for  a  transient  or  "standing  wave"  motion 
in  the  same  system  the  condition  for  a  shockless  motion  is  (see 
Mortell  and  Seymour  (1972a)) 

|MA(t)|  <  (i+a)  .  (4.11) 

where  A(t)  here  is  the  acceleration  level  defined  by  the 
initial  conditions.  The  formation  of  a  shock  is  determined  by 
the  induced  acceleration  level  in  the  gas  flow.  For  the  shockless 
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transient  motion  the  applied  and  induced  accelerations  have  the 
same  order  of  magnitude.  However,  for  the  resonant  forced 
motion,  the  induced  acceleration  has  the  same  order  of  magnitude 
as  the  square  root  of  the  applied  acceleration.  With  this 
observation  the  results  (4.10)  and  (4.11)  are  in  harmony.  In 
a  study  of  a  radiating  gas,  Eninger  (1971)  found  numerically  a 
critical  damping  which  prevents  shocks  at  resonance.  For  his 
analysis  damping  enters  as  a  weighted  integral  of  F  in  the 
equivalent  of  (2.27). 

Vi'hen  the  resonant  motion  is  shockless,  linear  theory  is  a 
uniformly  good  approximation  to  the  nonlinear  theory  provided 
the  piston  acceleration  is  sufficiently  small,  (see  figure  [4  ]). 
By  setting  v=0  in  (2.27),  the  linear  solution  is  given  by 


FL(n)  =  G(n)/y 


(4.12) 


and  hence,  by  (2. 27)  , 

FlnT"  ■  1 


^|F’(n)l 


(4.13) 


By  differentiating  (2.27)  and  setting  F"(ri)=0  we  find  that 
|F'(n)|  £  ^^(n)  where  A*^(n)  is  given  by  (3.24).  If  now  the 
applied  acceleration  is  small,  in  the  sense  that 


4vH *  (n) 


«  1  , 


(4.14) 


then 
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«  1  ,  by  (4.14). 


The  inequality  (4.9)  defines  a  critical  acceleration  level 
which  provides  a  sufficient  condition  on  the  applied  rate  to 
ensure  a  response  of  the  gas  which  is  continuous.  Numerical 
integration  of  the  equations  shows  that  shockless  solutions 
exist  for  piston  accelerations  greater  than  the  critical  one, 
so  that  inequality  (4.9)  is  conservative,  as  may  be  expected. 
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5.  Comparison  with  Experiments 

The  theoretical  predictions  of  the  analysis  presented  are 
compared  with  some  experimental  measurements  made  by  Sturtevant 
(1972).  His  set  up  consists  of  a  tube  of  length  132.5  inches 
with  an  inside  diameter  of  3.0  inches  which  contains  air  (y=1.4). 
At  one  end  of  the  tube  is  a  piston  which  is  displaced  sinusoidally 
with  an  amplitude,  normalized  against  the  length  of  the  tube,  of 
.0147  .  The  experiments  we  are  concerned  with  here  have  the 
two  configurations 

i)  the  far  end  of  the  tube  is  closed; 
ii)  the  far  end  of  the  tube  has  a  hole  in  it  whose  area 
is  small  compared  to  the  rea  of  the  end. 

For  case  i)  we  are  concerned  with  measurements  of  the  following 
quantities,  at  the  closed  end,  for  values  of  the  piston  fre¬ 
quency  around  the  fundamental: 

a)  the  absolute  maximum  and  minimum  of  a  normalized 
pressure  waveform, 

b)  the  pressure  immediately  before  and  after  the  shock 
jump. 

As  a  consequence  of  these  readings  the  values  of  the  frequency 
corresponding  to  the  lower  and  upper  ends  of  the  resonant  band 
are  available. 

Figure  [  5  ]  shows  the  comparison  between  inviscid  theory 
(i=0)  and  experimental  "response  curve"  of  Sturtevant.  In  this 
case,  our  theory  is  equivalent  to  that  in  section  4  of  Chester 
(1904).  In  computing  the  theoretical  curves  in  Figure  [6  ]  a 
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1  -k 

value  of  i  =  (i=.08)  is  chosen  so  that  the  shock  strength 

1+k 

exactly  at  resonance  C<S=0)  is  equal  to  the  observed  strength. 

The  theory  predicts  that,  for  6=0  ,  the  pressure  immediately 
before  the  shock  is  the  negative  of  the  pressure  immediately 
afterwards,  which  is  not  the  experimental  result.  Thus  the 
theoretical  and  experimental  curves  do  not  coincide  at  6=0  , 
even  with  our  choice  of  k  .  We  should  also  bear  in  mind  that 
for  the  conditions  of  the  experiment  the  small  rate  condition 
is  only  marginally  satisfied  (see  the  coimnent  at  the  end  of 
section  2) .  The  experiments  show  that  at  resonance  the  maximum 
pressure  exceeds  the  pressure  ahead  of  the  shock  and  the  pres¬ 
sure  behind  the  shock  exceeds  the  minimum  pressure.  This  is 
not  predicted  by  inviscid  theory,  but  is  a  property  of  the 
solution  of  the  equation  with  damping.  Another  point  to  note 
is  that  the  absolute  maximum  of  the  pressure  occurs  about  101 
to  the  right  of  the  resonant  frequency  while  the  absolute  min¬ 
imum  occurs  about  51  to  the  left.  For  the  experimental  conditions 
here  boundary  layer  damping  has  little  effect.  Nevertheless, 
if  the  other  damping  mech?''^sms  were  absent,  boundary  layer 
damping  could  prevent  a  if  the  radius  of  the  tube  and  the 

amplitude  of  the  piston  were  appropriately  adjusted. 

An  interesting  point  is  that  the  amount  of  damping  required 
to  get  the  shock  strength  correct  for  6=0  has  a  negligible 
effect  on  the  width  of  the  resonant  band.  This  might  seem  sur¬ 
prising  since  damping  decreases  the  shock  strength  which  in  turn 
determines  the  resonant  band.  The  result  can  be  understood, 
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qualitatively ,  by  considering  equation  {2.21)  and  bearing  in 
mind  the  definition  of  G(n)  given  by  (2.24).  When  there  is 
damping,  yfO  ,  the  system  defined  by  (2.27)  is  being  driven 
by  the  forcing  functior  G  whose  mean  is  non-zero  for  6f0  . 

The  increased  damping  is  then  counteracted  by  the  increased 
amplitude  of  the  effective  driver  G  . 

For  case  ii)  it  is  observed  that  for  particular  experimental 
conditions  there  is  a  critical  area  ratio  at  which  the  shock  in 
the  tube  disappears  for  all  frequencies.  If  we  interpret  the 
presence  of  a  small  hole  in  the  end  of  the  tube  as  a  means  of 
introducing  damping  into  the  system  then  the  prediction  of  the 
theory  agrees  qualitatively  with  experiment.  It  cannot  be  ex¬ 
pected  that  the  impedance  condition  (2.4),  as  introduced  in  the 
theory,  will  account  for  the  detailed  motion  of  the  gas  near 
the  orifice.  .Nevertheless,  it  seems  to  be  u.ssfal  in  predicting 
the  gross  features  of  the  motion. 

Curves  of  shock  strength,  S(i)  versus  impedance,  i  , 

were  plotted  for  various  values  of  the  piston  amplitude,  e  , 
Figure  [  7]  shows  there  is  a  linear  relation  between  S(0)-S(i) 
and  i  ,  which  is  independent  of  e  for  0£i<.2  when  .0l£e£.C2 
A  corresponding  plot  of  shock  strength  versus  area  ratio  would 
give  a  measure  of  the  effective  impedance  (or  effective  damping). 
The  linear  relationship  indicates  that  shock  strength  is  a  good 
measure  of  damping. 

This  result  can  be  understood  from  a  rough  analysis  of  the 
energy  balance.  When  i<<l  ,  the  result.s  of  section  3  indicate 
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that  if  the  amplitude  of  the  piston  is  e  ,  the  amplitude  of 

k  k 

the  response  is  O(e^)  ,  while  the  shock  strength  S  is  O(e^) 

The  balance  between  the  input  of  energy  due  to  the  piston  and 

the  loss  due  to  the  shock  and  radiation  from  the  end  is 

Ac =*  S^(i)  +  iBe  ,  (5.1) 

if 

where  A,  B  are  constants.  Since  S  =  0(e)  ,  (5.1)  can  be 
interpreted,  dividing  through  by  e  ,  as 

S(i)  -  S(0)  =  iB' 

where  S(0)  =  A*e^  ,  and  A',  B'  are  constants.  The  linear 
relationship  is  lost  when  shock  dissipation  is  no  longer  a 
major  effect.  As  the  piston  amplitude  (e)  decreases,  the 
point  at  which  the  curves  bifurcate  moves  towards  the  origin, 
so  that  the  linear  relation  holds  for  a  smaller  range  of  the 
impedance . 


We  sketch  a  regular  perturbation  procedure  for  the  deriva¬ 
tion  of  equation  (2.27)  directly  from  the  governing  equations 
(2.1) -(2. 5).  The  procedure  is  a  generalization  of  that  used  by 
Mortell  (1971b).  Ws  assume  an  expansion  of  the  form 

u 

u(t,x;e)  s  u^(t,x)  +  e  U2(t,x)  +  ... 
k 

e(t,x;e)  =  e  ej^(t,x)  +  e  +  ... 

0)  =  «»)(e)  “  w  +  0),  +  eWo  *  ••• 

o  X  d 

i  =  i(e)  »  i^  +  eig  +  . . .  , 

where  the  perturbation  parameter  is  the  amplitude  of  the  response. 

It 

On  noting  that  under  the  above  expansion  the  problem  at  O(e^) 
is  homogeneous,  we  find 

u^  »  ftn+WgX)  -  f(n'a)^x)  ;  -  =  f(n+a)^x)  +  f(n-WQx) 

where 

n  —  wt  ,  ^0  ^  2^  *  ^  ^  1,2,3,...  , 

2u) 

and  f(n)  is  an  arbitrary  function  ».lth  period  .  According 

k 

to  the  expansion,  is  0(1)  ,  so  that  e  n  <<  1  ,  which  is 

the  small  rate  condition. 

The  problem  at  0(e)  is  ncn-homogeneous ,  but  it  can  be 


integ-  .ted  to  give 
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2u2(t.x)  =  gj  *  ♦  2u’x(f;+f:)  *  ^  (fj-f=) 

n+w  X  n-w  X 

+  (Y+1)  |w^x(f_r+f^f;)  +  if2  (  £°ds  -  if;  j  £°dsj 

262(1, X)  =  gj  -  fg  ♦  2uix(f:-f;)  ♦  *  f^f.] 

n+w^x  n-w^x 

.  +  (y+1)£0^x(£_£;;^+£^£;)  +  I  £  ds  +  I  £  dsj 

where  the  subscripts  +,  -  indicate  the  argument  o£  the  function 
is  H+WqX  s  respectively,  prime  denotes  differentiation 

with  respect  to  the  argument,  and  gg  and  fg  are  arbitrary 
functions  associated  with  the  general  solution  of  the  homo¬ 
geneous  wave  equation.  The  boundary  conditions  at  this  order 
are 

Ug(t,0)  =  -i^e^(t,0)  and  Ug(t,l)  =  w^h'Cn)  . 

2(0 

Upon  using  the  conditions  that  £  and  gg  have  period  — °  , 
and  that  £  has  zero  mean  over  this  period,  these  boundary 
conditions  imply  that  £  satisfies 

Wo(Y+i)fCn)£’Cn)  +  2(o^£’(n)  +  2i^£(n)  =  Woh‘(n-(o^)  . 

If  we  now  note  that 

,,  _  1  „  Y'^l  -  M 

‘*'0  "  I  -  M 

and  define 

to'  =  w  to,  and  F  f  +  to./M  , 

1  o  1  i. 


the  equation  for  F  is 
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Zi.Ui  , 

nMFCn)F’Cn)  +  2i^F(n)  =  — ^  +  y  nh’ (n-w^)  .  (Al.l) 

This  agrees  with  equation  (2.27)  except  for  the  phase  of  h  . 
This  is  accounted  for  by  noting  that  the  parametrization  of  the 
S-wave  is  chosen  differently  here  (for  convenience  in  the  in¬ 
tegration)  than  in  the  body  of  the  text.  The  perturbation 
scheme  is  predicated  on  the  assumption  that  the  amplitude  of 
the  response  is  0(e)  .  There  is  then  an  implicit  assumption 

on  the  amount  of  damping  present,  and  hence  the  restriction 
k 

i=0(e  )  .  The  derivation  of  (2.27)  avoids  this  limitation  on 
its  range  of  applicability.  There  is  no  prior  assumption  on 
the  final  amplitude  and  hence  none  on  the  impedance,  i  .  Con¬ 
sequently,  from  (2.27)  it  is  seen  that  if  i=0  ,  nonlinearity 

k 

dominates  and  the  resulting  amplitudes  are  0(e)  ,  while  if 


i=0(l)  ,  i.e.,  damping  dominates,  then  the  amplitude  of  the 
response  is  0(c)  .  In  this  sense  equation  (2.27)  is  uniformly 
valid  in  the  damping  parameter,  whereas  (Al.l)  is  not. 
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Appendix  2. 

Here  we  prove  two  results  used  in  section  3: 

(i)  A  sufficient  condition  for  the  existence  of  a  continuous, 
periodic  solution  of  equation  (3.22),  for  all  6  ,  is  that 


>  max  [-4vH*Cn)] 

n 

(ii)  Given  pe(0,p^),  6e(6”,6'*')  and  F 

exists  an  n  e(n«,n  )  such  that 
s  &  o 

J  F(s)ds  =  ^  . 


(A2.1) 


defined  by  (3.29),  there 


(A2.2) 


Proof  of  (i) 

Using  the  notation  of  section  3  we  must  show  that  when 
(A2.1)  holds  the  separatices  and  Z"  pass  through  the  node 

,  i.e.  Z^(0^)  =  ^2(6^)  =  0  .  We  will  prove  the  result  for 
Z*  ;  the  argument  for  Z“  is  similar. 

Firstly,  since  the  curve  G(n)/p  is  the  isocline  Z'(n)=0 
and  Z(n)->-0'‘’  yields  the  isocline  2'(n)-*-+"  ,  for  6^<n<0j^  , 
the  separatrix  Z^  is  continuous  and  differentiable  in 
and  satisfies 

0  <  Z^(n)  1  max  [^;^]  . 

n  ^ 

In  particular  Z^(0^)  ^  0  .  We  show  that  3^(0^)  =  0  by  bound¬ 
ing  Z^(n)  above  by  a  function  Y(n)  which  has  the  properties 
Y(n)  >0  for  and  Y(0p  =0  .  Such  a  curve  bounds 
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dZ'*' 

above  if  <  Y'(n)  »  for  all  rie[6Q>0j_]  >  whenever  Z‘‘'=Y 

The  curve  Y(n)  =  has  these  properties  whenever  (A2.1) 

holds.  For  when  Z^’  =  Y  , 


dZ*' 

0 

3n  " 


^  n  IF  <  IHULni 
vY  2v  y 


which  holds  whenever  condition  (A2.1)  does.  Hence  since 
v(6j^)  =  0  ,  z^ce^)  =  0  . 

Proof  of  (ii) 

Defining 

y(n^)  =  J  PCs)ds  =  z^(s)ds  ♦  J  Z-(s)ds 


6 


and 


o  o  s 

we  wish  to  show  that,  for  a  given  6  ,  there  is  a  value  of 
ngECng.no)  such  that  y(ng)  =  g-  ,  where  0o<n2<0i<nQ<02  » 
ZqCDq)  =  ZgCHg)  =  0  .  We  first  note  that,  for  a  given  6  ,  y 
is  a  continuous  function  of  n  .  Since  is  continuous  in 

(00, Hq)  I  integration  of  (3,22)  yields 


n  n 

f  0  1  f  0 

z;!;(s)ds  =  i  G(s)ds  . 

J  °  F  J  f) 


Then  consider  y(n„)  = 


z;:;Cs)ds  +  z:(s)ds 
e  •’n. 


■o 

•e. 


iMGC=)ds  *  |  =(Z-(S) 


-  GCs) 


ds  . 


Gini 

-  U 


Thus,  since  Z  (n)  >  for  0^<n<6o  , 


y(n  )  >  -  ^ 

/  '■  'O-'  _  jj 


e. 


G  (s )  ds 


6 

S’  • 
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&  ! 

Similarly  yCHg)  1  5-  •  Therefore  by  the  continuity  of  y  , 
there  is  an  n„e(no»n-)  such  that 

S  O 

y(n3)  =  FCs)ds  =  «  . 

”0 
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Figure  1.  Construction  of  solution,  F  ,  using  integral  curves 
Shock  is  at  n 
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Figure  2.  F  =  Z  is  solution  at  edge  of  resonant  band,  5  =  .14 


,  0147 
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Figure  3.  F  =  Z  is  continuous  solution  outside  resonant  band. 


Figure  4.  Highly  danped  case;  linear  theory  a  good  approximation. 


Experimental  max  and  min  pressure 

Experimental  pressure  before  and 
after  shock 


Figure  6.  Theoretical  and  experimental  response  curves:  damped  case 


Shock  amplitude  for  i=0  minus  shock  amplitude 


